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ABSTRACT 

We attempt to characterize the observed variability of the magnetic helium-strong star a Ori E in 
terms of a recently developed rigidly rotating magnetosphere model. This model predicts the accumu¬ 
lation of circumstellar plasma in two co-rotating clouds, situated in magnetohydrostatic equilibrium 
at the intersection between magnetic and rotational equators. We find that the model can reproduce 
well the periodic modulations observed in the star’s light curve, Ha emission-line profile, and longitu¬ 
dinal held strength, confirming that it furnishes an essentially correct, quantitative description of the 
star’s magnetically controlled circumstellar environment. 

Subject headings: stars: individual (HD 37479) — stars: magnetic fields — stars: rotation — stars: 
chemically peculiar — stars: early-type — circumstellar matter 


1. INTRODUCTION 

As befits its status as the archetype of the helium- 
strong stars, characterized by elevated photospheric he¬ 
lium abundances, the B2Vpe star a Ori E is the most 
studied in its class. This is predominantly due to 
the rich phenomenology of variability manifested by 
the star, w hich exh i bits m odulations in its Ha Balmer 
emission (IWa.lhornl IToTl. its helium absorption-line 
strengths Pederseim fc_Thomsenlll97'i!h its photometric 
indices llHesser et al.ll 19761). its IJV continuum flux and 
resonance-line st rengths llSmith fe Groo tel200lJ) . its 6 cm 
radio emission jjLeonejfc Ilmanal 1 199. ‘ill . its linear po¬ 
larization lIKemn fc Hermanl ^177]) , and its longitudi- 
nal magnetic held strength llLandstreet fc Borrall978t) 
all varying on the same 1.19 dy period identihed 
with the stellar rotation cycle. Based on these exten¬ 
sive observations, an empirical picture has emerged fe.g. , 
iGroote fc Hungerl11982 : IBolton et al ] 5198 7: fShorej099^1 

of the star as an oblique-dipole magnetic rotator, whose 
held supports two circumstellar clouds situated at the 
intersections between magnetic and rotational equators. 


Following on from previous attempts at devising a 


theoretical basis for this picture (see, e.g., Nakaiima 

1985 

Shore fc Brown 1990; Bolton; 199J Short fc Bolton 

1994 

*), [Townsend & Owocki (2005. hereafter TO) have 


recently presented a new theoretical model for the distri¬ 
bution of circumstellar plasma around stars, like cr Ori E, 
that are characterized both by strong magnetic helds 
(such that the magnetic energy density is very much 
greater than the plasma energy density) and by signif¬ 
icant rotation. The model is built on the notion that 
plasma, flowing along hxed trajectories prescribed by 
rigid held lines, experiences an effective potential arising 
from a combination of the stellar gravitational held and 
the centrifugal force due to co-rotation. In the vicinity 
of points where the effective potential undergoes a local 
minimum, as sampled along a held line, the plasma can 
settle into a stable configuration of magnetohydrostatic 
equilibrium. 


These concep ts, _ originally advanc ed b y 

fMichel fc Sturrockl (I1974D and iNakaiimal (1985), 
are augmented in the TO treatment with a descrip¬ 
tion of the process Riling the wells established by 
local potential minima. In a manner similar to the 
magnetically confined wind shock (MC WS) model 
envisaged bv lBabel fc Montmerlei fiaazaH), radiatively 
driven wind streams from opposing footpoints collide 
at the top of closed magnetic loops, producing strong 
shocks. The post-shock plasma cools via radiative 
emission, and then either fal ls back to the star (see, 
e.g-.lud-Doula fe QwocOl2002l their Fig. 4), or — if the 
rotation is sufficient — settles into an effective potential 
well, where it continues to co-rotate in the circumstellar 
environment and accumulate over time. 

A strength of this rigidly rotating magnetosphere 
(RRM) model is its ability to predict quantitatively, 
via a semi-analytical approach, the relative circumstel¬ 
lar plasma distribution for arbitrary magnetic held con¬ 
figurations — not just for the simple axisymmetric case 
of a rotation-axis aligned dipole, to which magnetohy- 
drodynamical simulations have so far been restricted on 
grounds of computational tractability. In particular, ap¬ 
plication of the RRM formalism to an oblique-dipole 
model star leads to a specihc prediction (cf. TO) of 
a density distribution that is sharply peaked into two 
clouds, situated at the intersections between magnetic 
and rotational equators. Such a distribution coincides 
with the observationally inferred picture of cr Ori E; thus, 
it is clear that the RRM model holds promise for under¬ 
standing both this particular star, and related objects. 

This letter presents results from a preliminary investi¬ 
gation into whether the RRM model can simultaneously 
reproduce the spectroscopic, photometric and magnetic 
variability exhibited by cr Ori E. The following section 
discusses enhancements we have made to the model since 
the TO study; Section[3]then confronts the predictions of 
the model with the observed behaviour of the star. We 
discuss and summarize our hndings in Section @J 
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2. THEORETICAL DEVELOPMENTS 

Although the RRM treatment presented by TO is quite 
general, much of their analysis focuses on the specihc case 
of a spherical star, having a dipole magnetic held whose 
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TABLE 1 

RRM Model Parameters 


ft 

£* 

P 

i 

a 

(W) 


(deg) 

(deg) 

(R P ) 

0.5 

10- 3 

55 

75 

(-0.041,0.30,-0.029) 


origin coincides with the stellar origin. In the present 
analysis, we relax both of these restrictions. Specifically, 
we allow the star to assume a centrifugally distorted fig¬ 
ure, by defining its surface not as a sphere but as an iso¬ 
surface of the effective potential function d>(r) (cf. TO, 
their eqn. 11). This isosurface is set at <!> = 
where R p denotes the polar radius of the star, and 
here and throughout — the other symbols have the same 
meanings as in TO. 

We also allow for the possibility of a decentered dipole 
field, by introducing a vector a specifying the displace¬ 
ment of the magnetic origin from the stellar origin. Then, 
the position vector in the stellar reference frame, r, is re¬ 
lated to the magnetic-frame position vector r via 

r = Qf + a, (1) 

where the orthogonal matrix Q specifies a rotation by an 
angle f3 (the dipole obliquity) about the Cartesian y axis. 
In the TO treatment, a was set to zero, allowing simple 
expressions (cf. their eqn. 21) to be obtained relating 
the spherical polar coordinates in each reference frame. 
For non-zero a, such expressions do not exist, but the 
transformation between frames remains well defined. 

3. RRM MODEL 
3.1. General Considerations 

In Tabled we summarize the basic parameters for the 
RRM model of cr Ori E. The rotation rate Q = 0.5 fl c is 
an estimation based on the star’s 1.19 dy period, while 
the scale height parameter £* (cf. TO, eqn. 38) has 
been assigned a value typical to early-type stars. The 
choices of the obliquity, (3 = 55°, and observer inclina¬ 
tion, i = 75°, have been guided by certain characteristics 
of the spectroscopic and photometric observations pre¬ 
sented below. In particular, the unequal spacing of the 
light-curve minima of a Ori E, whereby the secondary 
eclipse follows the primary by ~ 0.4 of a rotation phase 
(see Fig. [TJ , restricts the possible RRM geometries to 
those satisfying the approximate relation f3 + i m 130°. 
A further indication that i > 75° comes from the fact 
that, at lower inclinations, the Ha variations — which 
take the form of a double S-wave (see Fig. 0) would 
exhibit a central spike when the two curves intersect; but 
no such spike is observed in the spectroscopy. 

The one remaining parameter of the model is the dipole 
offset vector a. The motivation for assuming a non-zero 
offset comes from the fact that the primary and sec¬ 
ondary eclipses have unequal depths; and that the two 
curves of the double S-wave are of unequal strength. The 
offset we select comprises (i) a displacement of the dipole 
center by 0.3 R p in a direction perpendicular to both 
magnetic and rotation axes, followed by (ii) a displace¬ 
ment of the dipole center by —0.05 R p along the magnetic 
axis. The first displacement produces a density contrast 



Fig. 1. - Observed (diamonds) and modeled (solid line) Strom- 
gren u-band light curves of a Ori E, phased on the star’s 1.19 dy 
rotation period. The error bars on the observational data are 
base d on conserva tive estimates of the scatter in the raw data 
fflesser et al.H1977l their Fig. 1). 

between the two clouds situated at the intersections of 
the equators, that allows us to achieve a good fit between 
the model and the observations. The second displace¬ 
ment has little effect on the circumstellar plasma distri¬ 
bution, but its introduction improves the fit to the mag¬ 
netic field observations. We note that IShort fc Boltonl 
lfT99l also suggested an offset dipole field, in order to 
explain asymmetries in the field observations. 

3.2. Photometry 

In Fig.Q we show the observed Stromgren u-band vari¬ 
ations of <7 Or i E, obtained via electronic extraction from 
Fig. 2 of lHesser et all ( 1977 1. Over these data we plot 
the synthetic light curve predicted by the RRM model. 
To calculate this curve, we assume that the photometric 
variations arise due to absorption of stellar flux by inter¬ 
vening material confined within the co-rotating magne¬ 
tosphere. Then, the observed u-band specific intensity 
X u , along a given sightline, is evaluated via a simple at¬ 
tenuation expression, 

I„=X*e _Tu , (2) 

where I* is the corresponding intensity at the stellar pho¬ 
tosphere, and t u is the u-band optical depth to the pho¬ 
tosphere. The variation of X* across t he stel lar disk is 
determined from a combination of the Ivon Zeipell 1192411 
gravity-darkening law (with the ansatz that the inten¬ 
sity varies_inproj3ortion to the bolometric flux) and 
the lEddingtonl ( 1926 1 limb-darkening law. We calcu¬ 
late t u from the density distribution predicted by the 
RRM model, under the assumption that the opacity in 
the magnetospheric material is spatially constant; this 
assumption corresponds, for instance, to bound-free ab¬ 
sorption processes in a medium having uniform compo¬ 
sition and excitation state. 

Integrating I u across the stellar disk, we arrive at the 
synthetic light curve shown in Fig. ^ Two free param¬ 
eters enter into this calculation. The normalization of 
I* is chosen so that the intrinsic (un-obscured) u-band 
brightness of the star is 6.69 mag. Likewise, the opacity 
and density are together scaled to achieve a match with 
the observed depth of the primary eclipse; this scaling 
corresponds to a maximal optical depth, over all sight¬ 
lines and rotation phases, of r„ = 1.63. 
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Fig. 2.— Observed (left) and modeled (right) time-series grayscales of the circumstellar Ho; emission of a Ori E, phased on the star’s 
1.19 dy rotation period. White indicates a 22% excess (in continuum units) over the background photospheric flux, and black indicates 
a 12% deficit; the grey levels are scal ed linearl y between the se extremes. The velocity axis is expressed in units of the star’s projected 
rotation velocity, usin i = 160kms —1 IGroote fc Hungeiil 19821) . 


The agreement between theory and observations is en¬ 
couraging; the synthetic light curve manages to cap¬ 
ture the timing, duration and strength of the eclipses. 
One obvious mismatch is that the RRM model is un¬ 
able to reproduce the emission feature at phase ~ 0.6, 
during which the star appears to brighten above its in¬ 
trinsic level. However, based on the fact that this feature 
is c onspicuously absen t from the u — b color variations 
(cf. iHesser et a tl ll977l) , we believe the feature may arise 
from some as-yet-unknown photospheric inhomogeneity, 
rather than from the circumstellar material. Further 
monitoring is required to test this hypothesis. 


3.3. Spectroscopy 

In Fig. |5J we plot the circumstellar H a variations ex¬ 
hibited by a Ori E, evaluated by subtracting a phase- 
dependent synthetic absorption profile from spectro¬ 
scopic observations of the star (25 echelle spectra, ob- 
tai ned during commis sio ning of the feros i nstrument; 
see lKaufer et al.ll999l and [Reiners et al.!ll20f)OH . The syn¬ 
thetic profiles are calculated usi ng an updated version of 
the model atmosphe re code o flHeberl rfl98! and a pho¬ 
tospheric model fsee lGroote fc Hungeitil997l) that takes 
into account the inhomogeneous helium surface distribu¬ 
tion. The resulting synthetic H(5 and Hy profiles are in 
very good agreement with the observations, on account 
of the absence of significant circumstellar contamination 
in these lines. Assuming that the synthetic Ha profile 
is also representative of the true (but unknown) pho¬ 
tospheric profile, we may interpret the observed-minus- 
synthetic difference profiles shown in Fig. |2] as an esti¬ 
mated circumstellar component of the total Ha flux. 

Alongside the observational data, we plot the corre¬ 
sponding predictions of the RRM model. We calculate 
the monochromatic specific intensity X\ at wavelength A, 
along a given sightline, via formal solution to the equa¬ 


tion of radiative transfer, 

I A =2*e- T * +S[l-e- T *]; (3) 

here, t\ is the monochromatic optical depth, and S is 
the source function, assumed to be constant and wave¬ 
length independent throughout the circumstellar envi¬ 
ronment. We model the variation of the photospheric 
intensity, across the stellar disk, using the same func¬ 
tion I* a dopted for the u-band synthetic photometry (cf. 
Sec. 13.211 : since we are ultimately interested only in the 
circumstellar component of the Ha profile, we do not in¬ 
clude any photospheric profile. For sightlines that do not 
intersect the disk, we set 2* to zero. 

To derive the optical depth t\, we assume an opac¬ 
ity proportional to the local plasma density; this choice 
reflects the fact that Ha absorption and emission, via 
radiative recombination, are density-squared processes. 
We further assume that the wavelength dependence of 
the opacity is a Gaussian of full width at half maximum 
AA, centered on the rest-frame wavelength, 6563 A, of 
the line. The modeled spectra are calculated by inte¬ 
grating I\ across the stellar disk and surrounding cir¬ 
cumstellar region. Three parameters are involved in the 
synthesis, which we tune by hand to achieve a good fit 
to the observations. This procedure leads to an intrin¬ 
sic line width A A = 30kms _1 , a source function S set at 
0.4 of the maximal photospheric intensity, and an opacity 
and density scaled such that the maximal optical depth 
achieved is t\ = 13.4. 

As with the photometry, there is good agreement be¬ 
tween theory and observations. The RRM model cor¬ 
rectly reproduces all of the qualitative features of the 
double S-wave Ha variations, including the asymmetry 
between the red and blue wings, the differing emission 
strengths at rotation phases ~ 0.25 and ~ 0.75, and the 
absorption features at phases ~ 0.0 and ~ 0.4. The phas¬ 
ing of the synthetic data is based on the photometric fit 
(cf. Fig. 0: thus, the small phase lag of the model, rela- 
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Fig. 3. — Observed (diamonds) and modeled (solid line) longi¬ 
tudinal magnetic field strength of a Ori E, phased on the star’s 
1.19 dy rotation pe r iod. The open diamonds come from the 
ILandstreet & B^ rrg[ J ML dataset, w hile the filled diamonds are 
taken from the Bohlend er et al.l <1987 ) dataset. 


tive to the observations, can wholly be att ributed to the 
uncer tainties in the rotation period (see lE.einers et alJ 
12000 ) compounded over the two decades separating the 
photometric and spectroscopic datasets. The contrast 
between the absorption features seen in the observations 
(strongest at phase ~ 0.4) and in the model (strongest 
at phase ~ 0.0) arises, most likely, due to recent changes 
i n th e star’s surface abundance distribution fsee lOrootei 
f2003 . his Fig. 3) that are not accounted for in our model. 

There are, of course, some discrepancies; for instance, 
the maximal separation of the S-wave curves is some¬ 
what smaller in the model (~ 2.75 v sin*) than in the ob¬ 
servations (~ 3.25 v sini). Nevertheless, given the many 
simplifying assumptions incorporated in the RRM model 
for a Ori E, it is quite striking how well it matches the 
observed spectroscopic variability of this star. 

3.4. Magnetic Field 

In Fig. 0 we show the time-varying longitudinal 
magnetic field strength o f a Ori E, as measured b y 
ILandstreet fc Borral ( 1978 1 and IBohlender et ~aT1 ( 1987 ): 
to maintain consistency with the spectroscopic data, 
we rephase the field d ata using the sam e period and 
ephemeris adopted hv lReiners et ah] (2 0001 ). Plotted over 
the observations is the field strength predicted by the 
RRM model. These s ynthetic data are calculated via 
an approach similar to iStibbsl ( 1950 ): the longitudinal 
(observer-directed) field is weighted by the local pho- 
tospheric intensity X* introduced previously, integrated 
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across the stellar disk, and then renormalized by the disk- 
integrated flux. The field geometry is specified by the pa¬ 
rameters given in Tabled and we choose a field strength 
such that the flux density is 11 kG at ±1 R p along the 
magnetic axis. Because the dipole is offset (cf. Sec. 13. 111 , 
the actual surface field strengths at the magnetic poles of 
the star are somewhat different than this nominal value. 

Once again, the agreement between theory and obser¬ 
vations is encouraging. The reduced chi-squared of the 
model, 1.79 for 20 degrees of freedom, is rather larger 
than t he value y 2 /n = 1.48 found bv IBohlender et all 
ltl987t) for a sinusoidal fit. However, these authors had 
the luxury of adjusting the phase of their fit, whereas in 
the present case the phasing is already constrained by the 
requirement that the primary light minimum be at phase 
0.0 (cf. Fig. We note that if the two outlier points 
at phases ~ 0.05 and ~ 0.15 are discarded, then the cor¬ 
respondence between model and observations improves 
greatly. Whether these points are indeed erroneous is an 
matter that should be resolved with further observations. 


4. DISCUSSION & SUMMARY 

In the preceding sections, we demonstrate how the 
photometric, spectroscopic and magnetic variability of 
a Ori E can be reproduced extremely well by a rigidly 
rotating magnetosphere model. This is persuasive evi¬ 
dence that the RRM paradigm furnishes an essentially 
correct, quantitative description of the star’s magneti¬ 
cally controlled circumstellar environment. 

Future studies can now focus on refining the model 
and its input parameters, beyond the exploratory treat¬ 
ment that we present here. Furthermore, the analysis 
can be extended to other observable quantities, such as 
radio emission, UV line strengths, and polarization. Ul¬ 
timately, based on our coarse exploration of parameter 
space, we expect that a fully optimized RRM model will 
be able to provide strong, independent constraints on the 
fundamental parameters (e.g., mass, radius) of a Ori E, 
allowing fr esh light to be she d on the historical uncer¬ 
tainty (see iHunger et al.1119891) concerning this excep¬ 
tional star’s distance and evolutionary status. 
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